Surface Scattering via Bulk Continuum States in the 3D Topological Insulator Bi 2 Se 3 
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We have performed scanning tunneling microscopy and differential tunneling conductance (dl/dV) 
mapping for the surface of the three dimensional topological insulator Bi2Se3. The fast Fourier 
transformation applied to the dl /dV image shows an electron interference pattern near Dirac node 
despite the general belief that the backscattering is well suppressed in the bulk energy gap region. 
The comparison of the present experimental result with theoretical surface and bulk band structures 
shows that the electron interference occurs through the scattering between the surface states near 
the Dirac node and the bulk continuum states. 

PACS numbers: 73.20.-r, 72.20.Dp, 72.25.-b 



Recent discovery of topological insulators has aroused 
great attention to the new state of quantum matter orig- 
inating from the surface state that forms a massless 
Dirac cone centered at time-reversal invariant momen- 
tum EHI1- Unlike the spin-degenerate Dirac fcrmions 
in graphcnc, this novel electronic state possesses heli- 
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FIG. 1: (color online) (a) STM and (b) differential conduc- 
tance images of Bi2Se3 surface at a sample bias voltage of 
—300 mV (10 nmx 10 nm). (c) Expanded image of the clover- 
shaped pattern of (b) (5 nmx 5 nm). (d) Fast Fourier transfor- 
mation image of (c). All data sets are taken at a temperature 
of 78K. 



cal spin textures protected by time-reversal symmetry, 
which creates the potential for spintronics devices with- 
out heat dissipation and for quantum computation ap- 
plications. The topological insulator phase has been pre- 
dicted to exist in a number of three-dimensional mate- 
rials, Bii.^Sb* 0, Bi 2 Se 3 , Bi 2 Te 3 , and Sb 2 Te 3 [H|. 
The spin-helical surface states have been experimentally 
observed by angle-resolved photoelectron spectroscopy 
(ARPES) 0,041, also with spin resolution I&EIEI- 

Among the binary chalcogenides, Bi 2 Se 3 is regarded as 
the most promising candidate owing to the single ideal 
Dirac cone residing in a wide band gap (0.3 eV). There- 
fore, significant effort has been applied towards the re- 
alization of spintronics devices using the quantum topo- 
logical transport in Bi 2 Se 3 [l4-17|. The magnetotrans- 
port measurements showed Shubnikov-de Haas oscilla- 
tions originating from the 3D Fermi surface, but no sur- 
face contribution to conduction has been observed in 3D 
topological insulators even at low carrier density [l5 - 17|. 

These observations pose the question of what limits 
the surface electron conduction in Bi 2 Se 3 . Band struc- 
ture calculations of Ref. predict the Dirac point of the 
surface state in Bi 2 Se 3 to be located below or close to the 
bulk valence band maximum, which may open a scatter- 
ing channel from the surface states to bulk continuum 
states 0|. In this Letter, we present scanning tunneling 
microscopy evidence of the near-surface scattering of the 
spin polarized surface states at the Dirac point in Bi 2 Se 3 
into spin-degenerate bulk continuum states. 

The scanning tunneling microscopy (STM) and spec- 
troscopy (STS) experiments were conducted at 78K and 
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FIG. 2: (color online) Experimental local density of states 
of Bi2Se3 surface acquired at 4.5 K (lower). Corresponding 
energy dispersion curve measured by ARPES for the surface 
state of Bi2Se3 is shown for comparison (upper). Both are 
shown in the same energy scale and their energy axes are 
fixed at the Dirac point. 



at 4.5K in an ultrahigh vacuum with a base pressure 
of lxlO -10 mbar using a low-temperature scanning tun- 
neling microscope (LT-STM, Omicron Nano Technology) . 
The STM images were acquired in the constant-current 
mode with a bias voltage V s applied to the sample. The 
dl/dV images were obtained simultaneously by record- 
ing the STM images using a lock- in technique. Bi 2 Sc 3 
samples were cleaved in situ at room temperature in 
an ultrahigh vacuum to minimize the surface contami- 
nation. For comparison, a photoemission spectrum was 
measured with a hemispherical photoclcctron analyzer 
(VG-SCIENTA R4000) at Hiroshima Synchrotron Radia- 
tion Center, Hiroshima University. The crystal structure 
of Bi 2 Se3 consists of weakly bonded 5-laycr (quintuple 
layer) blocks composed of Se-Bi-Se-Bi-Se atomic layers. 
Therefore, the topmost Se layer is exposed to the ultra- 
high vacuum. 

The STM topographic image and the dl/dV map ac- 
quired at 78K for the sample bias voltage V s = —300 mV 
are shown in Figs. 1(a) and 1(b). We find two differ- 
ent defect images with clover- and triangular-shaped pat- 
terns superimposed on the atomically resolved periodic 
structure spreading over the entire surface area in the to- 
pographic STM image. In the dl/dV image of Fig. 1(b) 
the clover-shaped pattern can still be seen, but not the 
triangular one. If we take a close look at the clover- 
shaped pattern in Fig. 1(c), we find the intensity modula- 
tion near the defect with a much longer wavelength than 
the surface lattice constant. A similar pattern has been 
observed in the pioneering STM study [l8( and identified 
as a defect state induced by an excess Bi atom substi- 



tuting the bottom Se site of the topmost quintuple layer. 
It was suggested that a resonance state is formed inside 
the bulk valence band due to the defect. Since the con- 
cept of topological insulator was not established at that 
time, we now need to revisit and reconsider the origin 
of the peculiar pattern in the STM image by taking into 
account the presence of the Dirac fermions with helical 
spin texture. 

In Fig. 1(d) we show the fast Fourier transformation 
(FFT) of the dl/dV image. Here the surface Brillouin 
zone (SBZ) boundaries are drawn from the obtained 
sharp FFT spots originated from the surface periodicity. 
The FFT image unambiguously tells us that an electron 
interference takes place as is seen from the flower shaped 
intensity distribution (Fig. 1(d)). The maximal intensity 
is seen to be along the TM direction of the scattering 
vector q, exhibiting a node in the TK direction. 

In order to further reveal the electron interference ef- 
fect around the defects in Bi2Se3, we have performed 
the experiment at a liquid He temperature (4.5 K). Fig- 
ure 2 shows the experimental local density of states 
(LDOS) measured as the differential tunneling conduc- 
tance dl/dV (STS) spectrum of Bi 2 Se 3 surface. We find 
a V-shaped STS spectrum with a minimum at —300 mV 
relative to the Fermi energy (Ep), which corresponds to 
the Dirac point (DP). The slope of STS spectrum is seen 
to become much steeper below —500 mV and above Ep- 
By comparing the present STS spectrum with the re- 
ported ARPES energy dispersion we find the bulk 
valence band (BVB) onset at —500 meV leading to the 
enhancement of LDOS below this energy. The shoulder 
at —100 meV coincides with the onset of bulk conduction 
band (BCB) as shown in Fig. 2. 

Figures 3(a)-3(i) represent the differential tunneling 
conductance dl/dV images from V s = —600 to —150 mV 
in the 50 nmx50 nm scanning area. Here the clover- 
shaped defect again appears with a brighter contrast with 
respect to the rest part of the images in a limited bias 
window between —450 and —300 mV (Figs. 3(c)-3(f)), 
which is close to and a little below the DP. However, it 
loses contrast at higher voltages (V s > —250 mV) and 
becomes even darker at V s < —500 mV (Figs. 3(a) and 
3(b)). 

To quantify the observed electron interference effect 
and determine the relevant q values, the FFT has been 
applied to the dl/dV images in Figs. 4(a)-4(i). At first, 
we observe a characteristic continuous intensity distribu- 
tion with sixfold symmetry inside the SBZ for V s = —450, 
—400, —350, and —300 mV, for which the clover-shaped 
patterns are observed in the corresponding dl /dV images 
(Figs. 4(c)-4(f)). Of the four FFT images, the intensity is 
the most prominent at V s = —400 and —350 mV, which 
coincides with the most distinct contrast in the clover- 
shaped pattern in Fig. 3(d) and 3(e). The characteris- 
tic six-fold pattern disappears above —250 mV, where 
only a weak intensity remains around the SBZ center 
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FIG. 3: (color online) (a)-(i) Differential conductance images 
of Bi2Se3 surface at 4.5K for several sample bias voltages from 
—600 mV to —150 mV in the area of 50 nmx50 nm. 



(Figs. 4(g)-4(i)). Below -450 mV, the FFT patterns be- 
come less anisotropic and their intensity is weakened as 
the bias voltage decreases. It follows from the FFT anal- 
ysis that the anisotropic interference is strongest between 
—450 and —300 mV, i.e., close to and slightly below the 
Dirac point. 

To gain deeper insight into the origin of the unexpected 
interference pattern, the cross-sectional FFT patterns 
traced along TM direction are shown in Fig. 4(j). For 
the whole range of sample bias voltages, we observe the 
largest intensity at g=0, which might originate from sev- 
eral extrinsic factors, such as a background noise whose 
modulation length would be much longer than the char- 
acteristic wavelength of the surface lattice periodicity and 
the electron interference. 

On the topological insulator surface, the spins arc 
locked with momenta, giving rise to helical Dirac 
fcrmions with lifted spin degeneracy. Owing to the 
strictly antiparallel spins at — k and k, surface electrons 
are protected from backscattering, and no standing wave 
pattern arising from the interference with the backscat- 
tered wave can be observed for a single Dirac cone. On 
the other hand, there might be other scattering channels 
within the surface state cone with a finite parallel spin 
component, which could contribute to the scattering with 
small q values of less than 0.2 A -1 (the maximal radius 
of the cone is 0.1 A -1 ). Although a detailed structure 
of surface state scattering is not resolved in the present 
data, the dominant contribution around q=0 apparently 



originates from the extrinsic effect because a finite inten- 
sity at q=0 persists and even enhances below —500 mV, 
where the surface states are merged into the bulk valence 
band. 

Let us now consider broad humps centered at q = 
±0.3 A -1 , which are especially pronounced in the sam- 
ple bias range from —450 to —300 mV as indicated by 
vertical bars in Fig. 4(j). The intensities in this cross 
sectional FFT profile are spread over a wide interval 
-0.5 < q < 0.5 A" 1 . Note that the FFT pattern is 
anisotropic, i.e., the intensity vanishes in the TK direc- 
tion, see Figs. 4(c)-4(f). The observed magnitude of the 
scattering vector, \q\ ~ 0.3 A -1 , is too large to be ex- 
plained by a scattering within the Dirac cone, especially 
near the DP. Besides, the profile is too broad to be ex- 
plained by a scattering between the surface states well 
localized in k|| . Also note that the intensity maxima 
of the hump are nearly independent of the bias voltage. 
Another important observation is that the hump almost 
disappears at 50 mV above the DP, which also excludes 
its origin from the surface state scattering. Since the 
scattering occurs at the near-surface defect, we need to 
consider both surface and bulk band structure to figure 
out a possible scattering channel other than the surface 
states of the Dirac cone. 

Wc have performed a first principles calculation of the 
surface states with a slab of seven formula units with 
a two-component self-consistent (within the local den- 
sity approximation) full-potential augmented plane wave 



method |22|. According to the theoretical band struc- 



ture the DP is energetically located slightly below the 
bulk valence band maximum in the middle of the TM in- 
terval, whereas around the T point the bulk band gap 
widens, so the states of the Dirac cone are energeti- 
cally well separated from the bulk states with small k|| 
vectors (Fig. 4(k)). The interference between a surface 
state with the lateral Bloch vector kjj and a bulk state 

with kjj caused by the scattering at the point defect re- 
sults in spatial oscillations of LDOS with the wave vector 
q = kjj — kjj . The possible s —> b transitions at the DP 

energy are seen to be centered at q w 0.3 A -1 in the f M 
direction. At the same time, in the TK direction such 
transitions arc forbidden in a wide energy interval: the 
valence band maximum is located much lower in energy 
than in the TM direction. This striking anisotropy ex- 
plains the negligible FFT intensity along TK line, which 
results in the flower-shaped distribution in Fig. 4. Since 
the bulk states are continuously distributed over the sur- 
face perpendicular momentum the interference pattern 
broadens, leading to a smeared hump in Fig. 4(j). Be- 
tween the valence band maximum and the conduction 
band minimum there exists a large portion of the upper 
cone that does not energetically overlap with bulk states. 
This is consistent with the abrupt disappearance of the 
hump above V s = —200 mV (Fig. 4(j)). This means that 
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FIG. 4: (color online) (a)-(i) Fast Fourier transformation of the dl/dV images in Figs. 3(a)-3(i) for respective bias voltages, 
(j) Cross-sectional profiles of the FFT images for (a)-(i) along TM . (k) Theoretical surface state dispersion (solid lines) and k|| 
projected bulk band structure (shaded area). An arrow shows possible surface to bulk scattering channels. 



the surface-bulk scattering is strongly suppressed above 
the DP, which supports our interpretation of the hump at 
q m 0.3 A -1 as due to the scattering into bulk states. A 
similar scattering may also occur via the bulk conduction 
band for V s > -100 meV with a smaller scattering vector 
\q\ < 0.2 A -1 , which is beyond the resolution limit of 
this measurement. Most probably, it is this process that 
gives rise to the triangular-shaped pattern in Fig.l, as it 
appears in the topographic image and disappears in the 
dl/dV map. 

In summary, we have shown by scanning tunneling mi- 
croscopy and differential tunneling conductance measure- 
ments that a strong surface state - bulk continuum scat- 
tering occurs near the Dirac node at the surface of the 
three dimensional topological insulator Bi2Se3, contrary 
to the belief that the surface state in the bulk energy 
gap around f point is protected from scattering ■ This 
may have implications for the construction of spintronic 
devices. In particular, it has been proposed to tune the 
DP at the top and the bottom surfaces of a thin film 
below and above the chemical potential in the dual-gate 
device configuration [23| ■ To suppress the scattering one 
would have to separate both the lower and the upper 
Dirac cone from the bulk states. Also, in order to realize 
the anomalous Hall effect by an appropriate magnetic 
doping, one must precisely tune the chemical potential 
right at the Dirac point to achieve the surface insulating 
state by opening the energy gap Q- Thus, we expect 
the present results to stimulate future band engineering 
in order to isolate the Dirac surface state from the bulk 
continuum and facilitate topological quantum transport. 
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